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If the surface tension of quark matter is low enough, quark matter is not self bound. At sufficiently 
low pressure and temperature, it will take the form of a crystal of positively charged strangelets in 
a neutralizing background of electrons. In this case there will exist, in addition to the usual family 
of strange stars, a family of low-mass large-radius objects analogous to white dwarfs, which we call 
"strangelet dwarfs". Using a generic parametrization of the equation of state of quark matter, we 
calculate the mass-radius relationship of these objects. 
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I. INTRODUCTION 

The matter that is directly observed in nature consists 
of atoms, whose nuclei are droplets of nuclear matter 
composed of up and down quarks. Nuclear matter is 
very stable: the most stable nuclei have lifetimes longer 
than the age of the universe. However, it has been hy- 
pothesized PH2] that nuclear matter may actually be 
metastable, and the true ground state of matter con- 
sists of a combination of roughly equal numbers of up, 
down, and strange quarks known as "strange matter". 
Strange matter is hypothesized to exist as (kilometer- 
sized) pieces, known as "strange stars" (reviewed in [3]), 
or as small nuggets, known as "strangelets" [3]. It has fur- 
ther been hypothesized that dark matter could be some 
form of quark matter, trapped in strangelets or strange 
stars before the era of nucleosynthesis [5HZ]- However, 
even if strange matter is not invoked as a dark matter 
candidate, there could still be a population of strange 
matter objects, from strange stars to strangelets, many 
of which would be relatively non-luminous. In this article 
we show that the masses and radii of such objects can ex- 
tend in to the range expected for planets. Recent surveys 
such as the Microlensing Observations in Astrophysics 
(MOA) and the Optical Gravitational Lensing Experi- 
ment (OGLE) to detect such low mass non-luminious 
low mass objects by gravitational lensing have yielded 
interesting results [5] . The hypothetical compact objects 
we predict could be detected by such methods and such 
surveys could place stringent bounds or perhaps hint at 
their possible existence. 

It is generally assumed that strange stars are compact 
objects, with sizes in the 10 kilometer range [4J, end- 
ing at a sharp surface of thickness ~ 1 Fermi, perhaps 
with a very thin electrostatically suspended nuclear mat- 
ter crust [&Mll|. However, if the surface tension a of 
the interface between quark matter and the vacuum is 
less than a critical value a CIit (of order a few MeV/fm 2 in 
typical models of quark matter) then large strangelets are 
unstable against fission into smaller ones [T21 - H4] . and the 
energetically preferred state is a crystal of strangelets: a 



mixed phase consisting of nuggets of positively-charged 
strange matter in a neutralizing background of electrons. 

In this "low surface tension" scenario, strange stars 
are not self-bound: they require gravitational attraction 
to bind the strangelets. Stars made of strange matter 
are then qualitatively similar to those made of nuclear 
matter: in each case the mass-radius relation has two 
branches, one compact and the other diffuse. For nu- 
clear matter, the compact branch contains neutron stars, 
which consist of gravitationally bound nuclear matter, 
with an outer crust that is a crystal of nuclei in a back- 
ground of electrons; the diffuse branch contains white 
dwarfs, which are a gravitationally bound cold plasma of 
nuclei (ions) and electrons, forming, at sufficiently low 
temperature, a crystalline structure (see, e.g., [15 ). For 
strange matter with a low surface tension, there are simi- 
larly two branches. The compact branch contains strange 
stars with a crust that consists of strangelets in a back- 
ground of electrons; this "strangelet crystal crust" was 
studied in Ref. [14] . In this paper we study the diffuse 
branch, which has no core of uniform quark matter: these 
stars consist entirely of strangelets in a background of de- 
generate electrons, so by analogy with white dwarfs we 
call them strangelet dwarfs. 

The strangelet-crystal phase is a charge-separated 
phase. Charge separation is favored by the internal en- 
ergy of the phases involved, because a neutral phase is al- 
ways at a maximum of the free energy with respect to the 
electrostatic potential (see [161 [Tf] ; for a pedagogical dis- 
cussion see [18). The domain structure is determined by 
competition between surface tension (which favors large 
domains) and electric field energy (which favors small 
domains). Debye screening plays a role in determining 
the domain structure, because it redistributes the elec- 
tric charge, concentrating it in the outer part of the quark 
matter domains and the inner part of the surrounding 
electron gas, and thereby modifying the internal energy 
and electrostatic energy contributions. Our parameteri- 
zation ([!} of the electrostatic properties of quark matter 
is generic, but is not appropriate for strangelets in the 
color-flavor locked (CFL) phase [19 , which is a degener- 
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ate case requiring separate treatment (see Sec. IV A I 



To obtain the M(R) relation of strangelet dwarfs, we 
solve the Tolman Oppenheimer Volkoff equation [2"0"]I21|. 
using the equation of state of the mixed phase. We obtain 
the equation of state by assuming that the strangelet 
lattice can be divided into unit cells ("Wigner-Seitz cells") 
and calculating the pressure of a cell as a function of its 
energy density. Our approach is similar to that used in 
previous studies of the strangelet crystal |131 [H] (except 
that in this paper we include electron mass effects) and 
in studies of mixed phases of quark matter and nuclear 
matter in the interior of neutron stars |22) . 

The main assumptions that we make are: 

1) We assume that the strangelets in the plasma form 
a regular lattice of Wigner-Seitz cells, which we treat as 
rotationally invariant (spherical) . In reality the cells will 
be unit cells of some regular lattice. We do not consider 
lower-dimensional structures (rods or slabs) because in 
Ref. [2] we found that such structures were never ener- 
getically favored. 

2) Within each Wigner-Seitz cell we use a Thomas-Fermi 
approach, solving the Poisson equation to obtain the 
charge distribution, energy density, and pressure. This 
is incorrect for very small strangelets, where the energy 
level structure of the quarks becomes important [23] [24] . 

3) We treat the interface between quark matter and the 
vacuum as a sharp interface which is characterized by a 
surface tension. We assume there is no charge localized 
on the surface. (Thus we neglect any surface charge that 
might arise from the reduction of the density of states of 
strange quarks at the surface [25H28] .) 

4) We neglect the curvature energy of a quark matter sur- 
face |29|, [30] , so we do not allow for "Swiss-cheese" mixed 
phases, in which the outer part of the Wigner-Seitz cell 
is filled with quark matter, with a cavity in the center, 
for which the curvature energy is crucial. 

5) We work at zero temperature. 

In our calculations we use units h = c = eo = 1, so 
a = e 2 /(47r) w 1/137. 



II. PHENOMENOLOGICAL DESCRIPTION OF 
QUARK MATTER 

We use the fact that in most phases of quark matter 
the chemical potential for negative electric charge fi e is 
much less than the chemical potential for quark number 
/i. This allows us to write down a model-independent 
parameterization of the quark matter equation of state, 
expanded in powers of fJ, e /fi |13| . 



(1) 



Note that the contribution of electrons to the pressure of 
quark matter is O(fif), and is neglected. This is a very 
good approximation for small strange quark mass, which 
corresponds to small riq . (For the largest value of Uq that 
we study, /i e in neutral quark matter is close to 100 MeV, 
and the assumption is still reasonable.) 



As noted in Sec. |T] we assume that the interface be- 
tween quark matter and vacuum has a surface tension cr, 
and we neglect any curvature energy. 

The quark density n and the electric charge density 
9qm ( m un its of the positron charge) are 
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So in uniform neutral quark matter the electron chem- 
ical potential is ^ eutral = Uq/Xq- Eq. |IJ) is a generic 
parametrization if ^n cutral <§; ^ which is typically the 
case in three-flavor quark matter. 

The bag constant enters in Po (/•*)) an d we will fix it by 
requiring that the first-order transition between neutral 
quark matter and the vacuum occur at quark chemical 
potential /x cr it, i- e - p(Mcrit, ^ GUtral ) = 0. Because we are 
assuming that the strange matter hypothesis is valid, we 
require /i cr it < 310 MeV, since at fi 310 MeV there is 
a transition from vacuum to neutral nuclear matter. In 
this article we will typically use /i cr it = 300 MeV. The 
value of /i inside our quark matter lumps will always be 
very close to /i cr it, so we can also expand in powers of 
A* — Merit, and write 



PQM^Me) ~ n (ji - Merit) + §X(M- Merit) 2 
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A quark matter equation of state can then be expressed 
in terms of 6 numbers: /i cr it, the charge density Uq and 
charge susceptibility Xq evaluated at fi = /i cr ; t , the quark 
number density n and susceptibility x evaluated at fi — 
/i cr it, and the surface tension a. 

We will restrict ourselves to values of the surface ten- 
sion that are below the critical value [T3J 



cr crit = 0.1325 — = 0.1325 
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where Ad is the Debye screening length in quark matter 
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If the surface tension is larger than cr cr it then the ener- 
getically favored structure at low pressure will not be a 
strangelet crystal, and there will be no strangelet dwarfs. 
Rough estimates of surface tension from the bag model 
are in the range 4 to 10 MeV/fm 2 [32 [35], and for typical 
models of quark matter, er cr it is of order 1 to 10 MeV/fm 2 
[13] . so it is reasonable to explore the possibility that 
strange quark matter could have a surface tension below 

Cent- 



A. Specific equations of state 

When we show numerical results we will need to vary 
Hq and Xq over a range of physically reasonable values. 
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To give a rough idea of what values are appropriate, 
we consider the example of non-interacting three-flavor 
quark matter, for which Uq and Xq become functions of 
fi and the strange quark mass m s , while po is in addition 
a function of the bag constant B. Expanding to lowest 
non-trivial order in m„ 



Po(p) 
n Q (n,m s ) 



9/x 4 
12^2 

2tt 2 
2n 2 
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We emphasize that these expressions are simply meant 
to give a rough idea of reasonable physical values for Uq 
and Xq ■ Our treatment does not depend on an expansion 
in powers of m s . To tune the transition between neutral 
quark matter and the vacuum so it occurs at /i = /i crit 
(see previous subsection), we set B so that po (Merit) = 

2 n 2 Q (Merit )/Xq Ocrit)- 

In the regions between lumps of strange matter, we 
will assume that there is a degenerate electron gas, whose 
pressure, and charge density in units of e, are 



Pe- (Me) 
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where fi^ = fcf. + m 2 . Note that at low pressures this 
is more accurate than the electron gas equation of state 
used in Ref . [TJ] , where the electron mass was set to zero. 



III. EQUATION OF STATE OF STRANGELET 
CRYSTAL 



A. Wigner-Seitz cell 

Following the approach of [11], we analyze a spherical 
Wigner-Seitz cell of radius i?ceih with a sphere of quark 
matter at the center of radius R. We use the Thomas- 
Fermi approximation to calculate /x e (r), 



V 2 /i e (r) = — 4:Traq(r) 



(8) 



where q(r) is the electric charge density in units of the 
positron charge e, and fi e is the electrostatic potential 
divided by e. 

The boundary conditions are that there is no electric 
field in the center of the cell (no 5-function charge there) , 
and no electric field at the edge of the cell (the cell is 
electrically neutral) , 
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We also need a matching condition at the edge of the 
quark matter. Since we assume that no charge is localized 
on the surface, we require continuity of \x e and its first 
derivative (the electric field) at r = R. 

The value of fi inside the strange matter will be slightly 
different from /i cr j t because the surface tension com- 
presses the droplet. To determine the value of /i, we 
require the pressure discontinuity across the surface of 
the strangelet to be balanced by the surface tension: 



PQM (M, Me (R) ) ~ Pe- (Me (R) ) 
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Once these equations are solved, we can obtain the 
equation of state of matter made of such cells. The total 
energy of a cell is 



E = 4tt 



r 2 dr {^n(n e ) - i/i e <j QM (^ e ) - _P qm (m, Me) 
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(11) 

The —^/j, e q terms in (111 come from combining — fi e q 
(from the relationship between energy density and pres- 
sure) with the electric field energy density +^i e q. The 
pressure of the cell is simply the pressure of the electrons 
at the edge of the cell, 



Pccll = Pe- (Me(-Rccll)) • 



The total number of quarks is 



N = Att r 2 dr n(/j,, \i e 



(12) 
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The volume of the cell is V = (4/3)7ri? co ii 3 . 

By varying R and i? C eii we generate a two-parameter 
family of strangelets. However, there is really only a 
single-parameter family of physical configurations, pa- 
rameterized by the external pressure p cc \\. On each line of 
constant p cc u in the (i?, i? C eii) parameter space, we must 
minimize the enthalpy per quark, 



h = 



N 



(14) 



to find the favored value of R and i? CG n. We assume zero 
temperature so h is also the Gibbs free energy per quark. 

We now have a well-defined way to obtain the equation 
of state of the mixed phase of quark matter, namely the 
energy density e — E /V as a function of the pressure 



B. Numerical solution 

Inside the quark matter, the solution to the Poisson 
equation Q that obeys the boundary condition at the 
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origin is 



1. Constant potential approximation 



A ■ W r \ 
, — smh( — ) 
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where A will be determined by matching conditions. 

In the degenerate electron gas region outside the 
strange matter, from ^ and (|8| the Poisson equation 
becomes 



V 2 /x e (r) 
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which must be solved numerically. For a given value of A 
we find from (15) the value and slope of /i e (r) at r = R, 
and use these as initial values to propagate )J. e (r) out to 
r = i? co n using ( |16[ ) . We vary A until we obtain a solution 
that obeys the boundary condition of no electric field at 
the edge of the cell. 



C. Low-pressure approximations 

If the pressure is not too high, the strangelet crys- 
tal consists of large Wigner-Seitz cells (R cc i\ 3> R)- In 
this regime one can obtain approximate analytic expres- 
sions for the equation of state of the crystal by assum- 
ing that the electrons have a roughly constant density 
outside the strangelet. We give these expressions below, 
and in later sections we use them to calculate mass ra- 
dius relations for large strangelet dwarf stars. However, 
we expect these approximations break down at ultra-low 
pressures, when the cell size becomes so large that screen- 
ing cannot be ignored, and the electrons are clumped 
around the strangelets, forming atoms, rather then being 
roughly uniformly distributed between the strangelets. 
This will happen when i? ce ii approaches the Bohr ra- 
dius ao = l/(am e ). i.e. when p ce u ^ a 5 Z 5 ^ 3 m^ « 
(10~ 12 MeV 4 )Z 5 / 3 . At these ultra-low pressures one 
should use an atomic matter equation of state: we do 
not do this, since we expect it will only affect a very 
small surface layer of the star, without any appreciable 
effect on the mass-radius relationship. 

The equation of state e(p cc u) is found by writing the 
energy density e of the cell and its pressure as a function 
of the size of the cell. For now we will treat the size R 
and charge Z of the central strangelet as unknowns; later 
we will estimate their values. 

Since the pressure inside a large cell is very low the en- 
ergy density of the quark matter is approximately n/z cr it, 



i? 3 

e w n/i cr it — — 3 



(17) 



To obtain the pressure at the edge of the cell we need to 
estimate the density distribution of the electrons outside 
the strangelet. 



The simplest approximation is to ignore screening, tak- 
ing the electron Fermi momentum kp e to be independent 
of r outside the strangelet (Sec. I of Ref. [33]) • Imposing 
neutrality of the cell fixes the Fermi momentum of the 
electrons, 
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Using ( 17 1, we obtain the equation of state e(p ce ii) of the 
strangelet crystal 



Merit n ~ 



4(fc Fe i?) 3 



9ttZ 



(19) 



where we use ([7]) to relate the electron Fermi momentum 

to Pcell. 

Because the constant potential approximation gives a 
fairly simple expression we can use it to understand how 
the strangelet crystal EoS depends on the parameters of 
the quark matter EoS, and hence how the M(R) curve 
for strangelet dwarf stars depends on those parameters. 
Note that in ( 19 1 the dependence of the energy density on 



the pressure is via a universal and monotonically increas- 
ing function kp e (p) ; dependence on the quark matter pa- 
rameters enters via the factor that multiplies this func- 
tion. To make the dependence on quark matter parame- 
ters explicit we use results for R and Z from Sec. |III C 3| 
below, and rewrite ( 19 ) for the EoS of the strangelet crys- 
tal as 



e(Pceu) 

s 



S (k Fe (p ccU )y , 

Merit n 

3TT 2 n Q £,(xo(a)) 



(20) 



where all dependence on the quark matter parameters 
comes through the prefactor S, which has units of en 
S can be explicitly obtained using ( 25 I 



ergy. 



(27 



and 

K Of S 



(28 I for the £ function. One could informally thin 
as a "softness" parameter of the strangelet crystal EoS: 
as S increases, the pressure becomes a more slowly-rising 
function of energy density. We expect that softer equa- 
tions of state will yield smaller stars with lower maximum 
masses. In Table |T] we give the value of S for a range of 
values of the parameters of the underlying quark matter 
EoS. 

At low enough pressures, the electrons become nonrel- 



kp e /(15ir 2 m e ), and (19) simplifies 



ativistic. Then p ce u 
to an analytic expression for the equation of state, 



4i? 3 /1257r o 

^z\-^ m * 
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3/5 

nMcritPc4u 



(21) 



This is a reasonable approximation when kp e < m e , i.e. 
when p C eii < m^/(247r 2 ) 0.0003 MeV 4 . However, as we 
will see below, at the very lowest pressures the constant 
potential approximation becomes inaccurate. 
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2. Coulomb potential approximation 

We can improve on the constant potential approxima- 
tion by including the Coulomb energy of the electrons in 
the calculation of the pressure. The equation of state is 
still given by ( 19 1 , but now the relationship between p cc \\ 



and kp e is modified by the addition of a Coulomb energy 
term (Ref. [33], (5)), yielding 



Pcell = Pe 



a / 
~ 5 \18ir 7 ) 



Z 2 \V3 ,, 4 



"Fe 



(22) 



Unlike the constant potential approximation, the 
Coulomb potential approximation gives an energy den- 
sity that goes to a non-zero value at zero pressure, 



ecoui(0) = Merit n 



2Z(am e Rf 
3^2 



(23) 



Comparing with (17 1 we see that this corresponds to the 
energy of cells with size of order l/(am e ) ~ 10~ 10 m. 
This is the energy density of a lattice of zero-pressure 
atomic matter with strangelets in place of nuclei, which 
is a reasonable guess for the low-pressure configuration of 
strangelets. We will therefore use the Coulomb approxi- 
mation as the low-pressure extension of our equation of 
state. As we will see, this leads to a "planet" branch 
in the mass-radius relation for configurations of strange 
matter. 



3. Radius and charge of strangelet at low pressure 

The low-pressure approximation expressions given 
above depend on the size R and charge Z of the strangelet 
at the center of a large cell. This is approximately an iso- 
lated strangelet, whose radius can be calculated by min- 
imizing the isolated strangelet free energy given in eqn 
(25) of Ref. H3J, 



Ag(x) = - 



3 x — tanh : 



35 

x 



(24) 



where x is the radius of the strangelet in units of A^> , and 



AiranQ\ 3 D 



(25) 



So the strangelet radius R as a function of the parameters 
of the quark matter equation of state is 



R = xqAd, where 



dAg 

dx 



(x ) = 



(26) 



We are interested in values of a up to 0.13, since for 
higher surface tension the strangelet crystal is no longer 
stable [13 . An approximate expression for the solution 
to (26 1, accurate to about 0.2% for a < 0.13, is 



approx 



15ct 



1/3 



2.174 a- 
1 - 3.982 5- 



(27) 



where the first term is the leading-order analytic expres- 
sion for xq in the limit of small a. 

The charge Z of the central strangelet is given by 
eqn. (17) of Ref. jT3], which can be written 



-irR 3 n Q aR/^D) 
-~ (x — tanh a;) , 



(28) 



where £ is a correction for the effects of screening inside 
the quark matter; it is an even function with £(0) = 1. 



IV. NUMERICAL RESULTS 

A. Range of parameters studied 

Our assumption that the strange matter hypothesis 
is valid requires that /Lt cr ;t must be less than the quark 
chemical potential of nuclear matter, about 310 MeV, 
so we fix /i cr it = 300 MeV. The value of /i inside our 
strange matter lumps will always be within a few MeV 
of /i cr jt i because if the surface tension is small enough to 
favor the strangelet crystal it will not cause significant 
compression. 

We will perform calculations for Ad = 4.82 fm and 
Ad = 6.82 fm, corresponding to Xq ~ 0.2/x^ rit (appro- 
priate for unpaired quark matter wol) and X 



0.1/A 2 



crit 

(appropriate for 2SC quark matter~[f3]). 

Typical values of Uq will be around 0.05^ cr itm^ 
and a reasonable range would correspond to varying m s 
over its physically plausible range, from about 100 to 300 
MeV. (To have strange matter in the star, m s must be 
less than /i cr it-) In this paper we use rig = 0.0445, 0.0791, 
and 0.124 fm -3 , which would correspond to m s = 150, 
200, and 250 MeV in (JsJ - 

There is another widely-discussed phase of quark mat- 
ter, the color-flavor locked (CFL) phase, but it is a de- 
generate case where rig = Xq = 0. CFL strangelets have 
a surface charge, but it does not arise from the mecha- 
nism studied here, Debye screening, and has a different 
dependence on the size of the strangelet [26 . We hope 
to study CFL strangelet matter in a separate work. 



B. Testing approximations to the equation of state 

In Fig. [I] we show the equation of state for the 
strangelet crystal, for critical quark chemical poten- 
tial /u cr it = 300 MeV, quark matter screening distance 
Ad = 6.82 fm, quark charge density parameter Hq — 
0.0791 fm -3 , and quark matter surface tension a = 
l.OMcVmV 2 . The dots were obtained numerically fol- 
lowing the procedure of Sec. |III B| The solid l ine is the 
Coulomb-potential approximation (Sec. Ill C 2 1 . On this 
plot the constant potential approximation (Sec. Ill C 1 1 
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TABLE I: Softness prefactor ( 20 1 of the strangelet crystal for various quark matter equation of state. The first two columns, Xd 
and Uq, specify the quark matter equation of state (|3jl (via (j5l). The third column gives the maximum surface tension for which 
a strangelet crystal will occur |4| . The last four columns givethe softness prefactor 5* for different values of the surface tension a 
(given in MeVfm -2 ) of the interface between quark matter and vacuum. 
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FIG. 1: Equation of state of the mixed phase (strangelet crys- 
tal) for strange matter with /x cr it ~ 300 MeV, \d = 6.82 fm, 
Uq = 0.0791 fm" 3 , a = l.OMeVfm" 2 . The dot s were ob- 
tained numerically following the procedure of Sec. |IIIB| Th e 



solid line is the Coulomb-potential approximation (Sec. Ill C 2 l 
The dashed line is the non-relativistic electron (ultra-low pres- 
sure) limit (21 1. Above p « 20000 MeV 4 , uniform quark matter 
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FIG. 2: Equation of state of the mixed phase for the same pa- 
rameters as in Fig.[l] zoomed in on the low pressure region, and 
with the energy density divided by p 0,6 . The dots were obtained 
numerically following the procedure of Sec. |III B"| The Coulomb- 
potential approximation (Sec. Ill C 2 1 is the most accurate, fol- 
lowed by the constant-potential approximation (Sec. |III C 1[ ), 
and then the non-relativistic electron approximation (21 1. 



line would be indistinguishable from the Coulomb- 
potential line, so we do not show it. The dot-dashed line 
is the non-relativistic electron (ultra-low pressure) limit 
(21 1 of the constant potential approximation. Above 



p w 20000 MeV , uniform quark matter becomes favored 
over the mixed phase. On this very expanded logarith- 
mic scale, the Coulomb approximation appears reason- 
ably accurate up to pressures of order 1 MeV. 

To achieve more discrimination between the different 
approximations, we show in Fig. [2] a magnified version 
of the low-pressure end of the plot in Fig. [T] where we 
have divided out the non-relativistic scaling of the en- 
ergy density, e ~ p 3 ^ 5 . We can see that, down to the 
lowest pressures for which we can perform the numeri- 
cal Wigner-Seitz calculation of the equation of state, the 



Coulomb approximation gives the most accurate semi- 
analytic approximation, although the constant potential 
approximation is accurate to within about 10%. 

We then have to decide which approximation to use 
for lower pressures, where numerical calculations are not 
available. In the low-pressure limit, the Coulomb approx- 
imation to e(p) tends to a fixed value, while the constant 
and nonrelativistic approximations to e(p) tend to zero 
as p 3 / 5 . So in Fig. ^ the Coulomb approximation will 
diverge at p -C 10~ 5 MeV 4 , while the constant and non- 
relativistic approximations will tend to the same constant 
value. As discussed in Sec. |III C 2| it seems reasonable to 
expect that at the lowest pressures there will be a crystal 
of "strange atoms", each consisting of electrons bound to 
a strangelet, and the Coulomb approximation gives a rea- 
sonable estimate of the energy density of such matter, so 
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FIG. 3: The full mass-radius curve for stars made of quark 
matter with the equation of state plotted in Fig. [T] using the 
Coulomb approximation | |22[ | to extrapolate to lower pressures. 
The compact branch contains strange stars with a strangelet 
crystal crust. The diffuse branch contains stars consisting en- 
tirely of strangelet crystal matter. Solid lines represent con- 
figurations that are stable; stability of the other branches is 
discussed in the text. 



at low pressure we will use the Coulomb approximation. 



C. Mass- radius relation of strange stars 

In Fig. [3] we show the full mass-radius curve for stars 
made of quark matter with the equation of state plotted 
in Fig. [T] The compact branch contains strange stars 
with a strangelet crystal crust. The diffuse branch con- 
tains stars consisting entirely of strangelet crystal mat- 
ter. It includes two segments: the lighter one is planets 
of dilute strange matter whose the mass increases with 
radius. This joins to the strangelet dwarf branch where 
the mass decreases with radius as the strangelet crystal 
is compressed by the pressure due to gravity. We use 
the numerically calculated equation of state (Sec 



III) 



except that at very low pressure (the planetary branch) 
the Wigner-Seitz cells become so large that our numer- 
ical methods break down, so as discussed in Sec. |III C 2| 
we use the Coulomb approximation ( 22 1 to extrapolate 
down to zero pressure. 

Fig. [3] shows the whole M(R) curve, not all of which 
corresponds to stable configurations. The usual stability 
criterion for stars |34| is that one radial mode becomes 
either stable or unstable at each extremum in the M(R) 
function. A stable mode becomes unstable at each ex- 
tremum where the curve bends counterclockwise as the 
central density increases; a stable mode becomes un- 
stable at each extremum where the curve bends clock- 
wise as the central density increases. However, Glen- 
denning et.al. [35 report that at some extrema there is 
no change in stability: the squared frequency of one of 
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are for a larger surface tension, a = 3MeVfm -2 . The 
dots were obtained using the full numerical equation of state 
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imation (Sec. |III C 2[ ), and the dashed lines use the constant- 
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FIG. 5: Mass-radius relation for strangelet dwarfs made of 
strangelet crystal matter, comparing different approximations 
to the equation of state. 



the fundamental radial modes may touch zero, but not 
change sign. We defer a detailed study of the stability of 
radial modes of strange stars to future work, and in Fig. [3] 
we show as "stable" (solid curves) the parts of the M(R) 
curve that both Ref. [31] and Ref. [35] agree are stable. 
We note that Ref. [35] is a study of stars that have a core 
of uniform strange matter surrounded by a crust of nu- 
clear matter: these are similar to the configurations along 
the dashed part of the mass-radius curve in Fig. [3] where 
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we have a core of uniform strange matter surrounded by 
a crust of strangelets, with a density discontinuity at the 
boundary. If Ref. [35] 's stability argument is correct and 
applicable to our stars, then some of these configurations 
may also be stable. In the remainder of this paper we 
will focus on the strangelet dwarf branch, which consists 
of a simple crystal of stranglets with no uniform core, so 
there is no controversy about the appropriate stability 
criterion. 



D. Mass-radius relation of strangelet dwarfs 

To investigate the sensitivity of the masses and radii of 
strangelet dwarfs to the parameters of the quark matter 
equation of state, we show in Fig.[4]and [5jthe strangelet 
dwarf part of the mass-radius curve, excluding the com- 
pact and planetary branches, for various values of the 
quark matter parameters. 

In Fig. [3] we explore the effects of varying the sur- 
face tension, and we compare the different approxima- 
tions to the equation of state. The upper curves are 
for the same equation of state as was shown in Fig. [l] 
and [2j the lower curves use a larger surface tension, 
a = 3MeVfm -2 . In both cases the solid curves are ob- 
tained from the Coulomb-potential approximation to the 
equation of state, and the dashed lines are obtained from 
the constant-potential approximation. The dots use the 
equation of state that is obtained numerically following 
the procedure of Sec. |III B| except that at very low pres- 
sures, where the numerical calculation becomes too diffi- 
cult, the Coulomb approximation is used. 

We see that, as one might have expected from Fig. [T] 
using the Coulomb approximation over the entire pres- 
sure range of the mixed phase yields reasonably accurate 
results. However, as noted in Sec. |III C 1| the constant 
potential approximation is still useful for gaining an un- 
derstanding of how the M (R) curve for strangelet dwarfs 
depends on the parameters of the EoS, because in the 
range of pressures that is important for strangelet dwarfs 
it gives a good indication of the M(R) curve. (At ultra- 
low pressures, relevant for the strange planet branch, this 
is no longer the case: one has to use the Coulomb approx- 
imation instead.) As discussed in Sec. Ill CI the con- 
stant potential approximation to the EoS can be written 
in terms of a "softness prefactor" S ( 20 1 . To understand 
how the M{R) curve in Fig. [4] changes with a, note that 
Xo(a) is a monotonically increasing function and £,(xq) 
is a monotonically decreasing function, so as the surface 
tension a increases at fixed values of the other parame- 
ters, the softness prefactor S of the strangelet crystal EoS 
increases (one can see this in Table |l| . Since the EoS is 
becoming softer, the M(R) curve moves down and to the 
left, giving smaller stars with a lower maximum mass. 

In Fig. [5] we explore the effects of varying the charge 
density parameter Uq in (|3l while keeping the other pa- 
rameters constant. As inTig. [4] solid lines are for the 
Coulomb approximation to the equation of state, dots 



are for the numerically calculated equation of state using 
the Coulomb approximation to extrapolate to the lowest 
pressures. We see that increasing rig yields heavier, larger 
strangelet dwarf stars. Again, this can be understood in 
terms of the consta nt p otential approximation and its 
softness prefactor 5* ( 20 1 . As riq increases, it causes S to 
decrease through two effects. Firstly via the explicit fac- 
tor of Uq in the denominator of (|20| , and secondly via the 
relationship p5| between a analr. The sensitivity of S 
to changes in rig can be seen in Table |lj for the two values 
of Uq used in Fig. |5]the values of S are near the extremes 
of its range in the parameter set we studied: S « 345 
and S fa 94 for n Q = 0.0445 and Vq = 0.124 respectively. 
Consequently, the M(R) curve for rig = 0.0445 is char- 
acteristic of a soft equation of state, with low radius at a 
given mass and a low maximum mass, whereas the M(R) 
curve for Uq = 0.124 is characteristic of a hard equation 
of state, with large radius at a given mass and a high 
maximum mass. 



V. DISCUSSION 

We have shown that, if the strange matter hypothesis is 
correct and the surface tension of the interface between 
strange matter and the vacuum is less than a critical 
value Q , there is at least one additional stable branch in 
the mass-radius relation for strange stars, corresponding 
to large diffuse objects that we call "strangelet dwarfs", 
consisting of a crystal of strangelets in a sea of electrons. 
This is easily understood, since if a < cr cr i t then uniform 
strange matter is unstable at zero pressure, and under- 
goes charge separation to a crystal of positively-charged 
strangelets surrounded by electrons, just as normal mat- 
ter at zero pressure is a mixed phase consisting of droplets 
of nuclear matter surrounded by electrons. Strangelet 
dwarfs are then the strange matter equivalent of white 
dwarfs. 

We emphasize that in this low-surface-tension scenario, 
strange matter is not self bound. Like nuclear matter, it 
is only bound by gravitational forces. Every strange star 
will have a strangelet crystal crust, and strangelet dwarfs 
are those strange stars that are "all crust". 

The natural production mechanism by which 
strangelet dwarfs might be produced is a collision 
between a strange star and another compact object. In 
such collisions, up to 0.03 Mq may be ejected (33], which 
is in the mass range we are predicting for strangelet 
dwarfs. There are two ways a collision could produce 
strangelet dwarfs. Firstly, part of the crust of the 
strange star might be ejected to become a isolated 
object, which would be a strangelet dwarf. Secondly, 
if a sufficiently light piece of the uniform quark matter 
core were ejected in the collision, it would be unable 
to exist on the compact branch, and would evaporate 
into a configuration on the diffuse branch. For example, 
for the equation of state studied in Fig. [3] the lightest 
compact configuration of strange matter is 0.0055 Mq. 
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A lighter piece of strange matter could only exist on the 
diffuse branch, and would spontaneously evaporate to 
become a strangelet dwarf. Strangelet dwarfs produced 
by these mechanisms could then bind gravitationally, to 
form heavier strangelet dwarfs. 

It should be noted that our proposed mechanism for 
the production of strangelet dwarfs is also a mecha- 
nism for creating a diffuse cosmic flux of strangelets 
("strangelet pollution") , which might be expected to con- 
vert all neutron stars to strange stars [37 . Although ob- 
servations of glitches and magnetar oscillations [35] seem 
consistent with some compact stars having nuclear mat- 
ter crusts, there remains some uncertainty. Crystalline 
phases of quark matter could allow strange stars to glitch 
|39| . and in our low-surface-tension scenario strange stars 
have crusts that could be hundreds of meters thick [2] . 
A cosmic flux of strangelets may seem unlikely but until 
it is ruled out experimentally (as may happen soon from 
the AMS experiment |JD]) it remains useful to analyze 
the full observational consequences of the strange matter 
hypothesis. 

Our analysis assumes that at any given pressure the 
strangelet crystal consists of the most energetically fa- 
vorable strangelet configuration (in terms of strangelet 
size and charge and cell size). However, other configu- 
rations will in general be metastable with long lifetimes. 
If one compresses a piece of strangelet crystal then the 
charge of the strangelets can readily change via absorp- 
tion or emission of electrons, but it is very difficult for 
the quark matter to rearrange itself in to strangelets of 
the now-energetically-favored size: it is more likely that 
the strangelets will stay the same size and the radial den- 
sity profile of the electrons will change. The sizes of the 
strangelets will be determined more by the history of the 
object than by the pressure. Taking this point further, it 
is quite possible to have a crystal consisting of a mixture 
of strangelets and ordinary nuclei, held apart by their 
electrostatic repulsion but also bound together in to a 
crystal by the degenerate electron gas that neutralizes 
them, forming a hybrid strangelet / white dwarf star. 

Detection of strangelet dwarfs requires an observation 
method that can find non-luminous objects with typical 



masses of 10 -5 to 10" 1 M Q and radii in the range 500 to 
5000 km. An example is gravitational microlensing sur- 
veys, such as those conducted by the Microlensing Obser- 
vations in Astrophysics (MOA) and the Optical Gravita- 
tional Lensing Experiments (OGLE) groups, which look 
for lensing events in the galactic bulge, and are capable of 
detecting Jupiter-mass objects. It is intriguing that such 
surveys now report the existence of an abundant pop- 
ulation of unbound distant planetary masses, suggest- 
ing that such objects may be twice as common as main 
sequence stars [8]. Although models of planet forma- 
tion indicate that mechanisms exist for unbinding planets 
through disk instabilities and planet interactions [H] , we 
suggest that a possible alternative is formation of strange 
dwarfs from matter ejected in strange star mergers. One 
would expect that sometimes a strangelet dwarf produced 
in a merger might be unable to escape the gravitational 
field of the remaining compact object, and this would 
explain the presence of dense planet-mass objects in the 
vicinity of compact stars. An example is the millisecond 
pulsar PSR J1719-1438, which has a Jupiter-mass com- 
panion whose inferred central density (p > 23 gem" 3 ) is 
far in excess of what is expected in a planet 02J. We 
expect that in the near future further light will be cast 
on this question, as microlensing surveys help us better 
understand the distribution of planetary mass compact 
objects and as strategies are devised to provide informa- 
tion about both mass and radius. 



Acknowledgments 

We thank E. Agol and F. Weber for useful discussions. 
This research was supported in part by the Offices of 
Nuclear Physics and High Energy Physics of the Office 
of Science of the U.S. Department of Energy under con- 
tracts #DE-FG02-91ER40628, #DE-FG02-05ER41375, 
and by the DoE Topical Collaboration "Neutrinos and 
Nucleosynthesis in Hot and Dense Matter", contract ^de- 
sc0004955. 



[1] A. R. Bodmer, Collapsed nuclei, Phys. Rev. D4 (1971) 
1601-1606. 

[2] E. Witten, Cosmic Separation of Phases, Phys. Rev. D30 

(1984) 272-285. 
[3] E. Farhi and R. L. Jaffe, Strange Matter, Phys. Rev. D30 

(1984) 2379. 

[4] F. Weber, Strange quark matter and compact stars, Prog. 

Part. Nucl. Phys. 54 (2005) 193-288, 

[astro-ph/ 0407155]. 
[5] A. R. Zhitnitsky, 'Nonbaryonic' dark matter as baryonic 

color superconductor, JCAP 0310 (2003) 010, 

[hep-ph/0202161]. 
[6] S. Banerjee et. al, Massive compact halo objects from the 



relics of the cosmic quark-hadron transition, Mon. Not. 
Roy. Astron. Soc. 340 (2003) 284, [astro-ph/0211560]. 

[7] M. M. Forbes, K. Lawson, and A. R. Zhitnitsky, The 
electrosphere of macroscopic 'quark nuclei': A source for 
diffuse MeV emissions from dark matter, 
arXiv:0910.4541. 

[8] T. Sumi, K. Kamiya, D. P. Bennett, I. A. Bond, F. Abe, 
C. S. Botzler, A. Fukui, K. Furusawa, J. B. Hearnshaw, 
Y. Itow, P. M. Kilmartin, A. Korpela, W. Lin, C. H. Ling, 
K. Masuda, Y. Matsubara, N. Miyake, M. Motomura, 
Y. Muraki, M. Nagaya, S. Nakamura, K. Ohnishi, 
T. Okumura, Y. C. Perrott, N. Rattenbury, T. Saito, 
T. Sako, D. J. Sullivan, W. L. Sweatman, P. J. Tristram, 



10 



A. Udalski, M. K. Szymanski, M. Kubiak, G. Pietrzyfiski, 

R. Poleski, I. Soszynski, L. Wyrzykowski, K. Ulaczyk, and [27 

Microlcnsing Observations in Astrophysics (MOA) 

Collaboration, Unbound or distant planetary mass 

population detected by gravitational microlensing, Nature [28 

(London) 473 (May, 2011) 349-352, [arXiv : 1105 . 3544]. 

C. Alcock, E. Farhi, and A. Olinto, Strange stars, 

Astrophys. J. 310 (1986) 261 272. [29 

M. Stejner and J. Madsen, Gaps below strange star crusts, 

Phys. Rev. D72 (2005) 123005, [astro-ph/0512144]. 

V. V. Usov, Low-mass normal-matter atmospheres of 

strange stars and their radiation, Astrophys. J. 481 [30 

(1997) L107 L110, [astro-ph/9703037]. 

P. Jaikumar, S. Reddy, and A. W. Steiner, The strange 

star surface: A crust with nuggets, Phys. Rev. Lett. 96 [31 

(2006) 041101, [nucl-th/0507055]. 

M. G. Alford, K. Rajagopal, S. Reddy, and A. W. Steiner, [32 

The stability of strange star crusts and strangelets, Phys. 

Rev. D73 (2006) 114016, [hep-ph/0604134]. [33 

M. G. Alford and D. A. Eby, Thickness of the 

strangelet-crystal crust of a strange star, Phys. Rev. C78 

(2008) 045802, [arXiv: 0808. 0671]. [34 

A. Y. Potekhin and G. Chabricr, Equation of state of 

fully ionized electron-ion plasmas. II. extension to 

relativistic densities and to the solid phase, Phys. Rev. E 

62 (Dec, 2000) 8554-8563, [astro-ph/0009261]. [35 

D. G. Ravenhall, C. J. Pethick, and J. R. Wilson, 
Structure of matter below nuclear saturation density, 

Phys. Rev. Lett. 50 (1983) 2066-2069. [36 

N. K. Glendenning, First order phase transitions with 

more than one conserved charge: Consequences for 

neutron stars, Phys. Rev. D46 (1992) 1274 1287. 

M. Alford, C. Kouvaris, and K. Rajagopal, Evaluating the [37 

gapless color-flavor locked phase, Phys. Rev. D71 (2005) 

054009, [hep-ph/0406137]. 

M. G. Alford, K. Rajagopal, and F. Wilczek, Color-flavor [38 
locking and chiral symmetry breaking in high 
QCD, Nucl. Phys. B537 (1999) 443 458, 
[hep-ph/9804403]. [39 
R. C. Tolman, Static solutions of Einstein's field equations 
for spheres of fluid, Phys. Rev. 55 (1939) 364-373. 
J. R. Oppenheimer and G. M. Volkoff, On Massive [40 
neutron cores, Phys. Rev. 55 (1939) 374-381. 
T. Maruyama, S. Chiba, H.-J. Schulze, and T. Tatsumi, 
Hadron-quark mixed phase in hyperon stars, Phys. Rev. [41 
D76 (2007) 123015, [arXiv: 0708. 3277]. 
J. Madsen, Shell model versus liquid drop model for 
strangelets, Phys. Rev. D50 (1994) 3328-3331, 
[hep-ph/9407314]. [42 
P. Amore, M. C. Birse, J. A. McGovern, and N. R. Walet, 
Colour superconductivity in finite systems, Phys. Rev. 
D65 (2002) 074005, [hep-ph/0110267]. 
J. Madsen, Intermediate mass strangelets are positively 
charged, Phys. Rev. Lett. 85 (2000) 4687 4690, 
[hep-ph/0008217]. 

J. Madsen, Color-flavor locked strangelets, Phys. Rev. 



Lett. 87 (2001) 172003, [hep-ph/0108036]. 

J. Madsen, Universal charge-radius relation for subatomic 
and astrophysical compact objects, Phys. Rev. Lett. 100 

(2008) 151102, [arXiv: 0804. 2140]. 

M. Oertel and M. Urban, Surface effects in color 

superconducting strangelets and strange stars, Phys. Rev. 

D77 (2008) 074015, [arXiv: 0801. 2313]. 

M. B. Christiansen and N. K. Glendenning, Finite size 

effects and the mixed quark-hadron phase in neutron 

stars, Phys. Rev. C56 (1997) 2858-2864, 

[astro-ph/9706056]. 

M. B. Christiansen and J. Madsen, Inhomogeneity scale 
from the cosmological quark-hadron transition, J. Phys. 
G23 (1997) 2039 2046. 

M. S. Berger and R. L. Jaffe, Radioactivity in strange 

quark matter, Phys. Rev. C35 (1987) 213-225. 

M. S. Berger and R. L. Jaffe, Erratum: Radioactivity in 

strange quark matter, Phys. Rev. C 44 (Jul, 1991) 566. 

E. E. Salpeter, Energy and Pressure of a 

Zero- Temperature Plasma, Astrophys. J. 134 (1961) 

669-682. 

J. M. Bardeen, K. S. Thorne, and D. W. Meltzer, A 
catalogue of methods for studying the normal modes of 
radial pulsation of general-relativistic stellar models, 
Astrophys. J. 145 (Aug., 1966) 505-+. 
N. Glendenning, C. Kettner, and F. Weber, A Possible 
new class of dense white dwarfs, Phys. Rev. Lett. 74 
(1995) 3519-3521. 

A. Bauswein, H.-T. Janka, R. Oechslin, G. Pagliara, 
I. Sagert, et. at, Mass Ejection by Strange Star Mergers 
and Observational Implications, Phys. Rev. Lett. 103 

(2009) 011101, [arXiv: 0812. 4248]. 

J. Friedman and R. Caldwell, Evidence against a strange 
ground state for baryons, Phys. Lett. B264 (1991) 
143-148. 

A. L. Watts and S. Reddy, Magnetar oscillations pose 

challenges for strange stars, Mon. Not. Roy. Astron. Soc. 

379 (2007) L63, [astro-ph/0609364]. 

M. Mannarelli, K. Rajagopal, and R. Sharma, The 

Rigidity of crystalline color superconducting quark matter, 

Phys.Rev. D76 (2007) 074026, [hep-ph/0702021]. 

J. Sandweiss, Overview of strangelet searches and Alpha 

Magnetic Spectrometer: When will we stop searching?, 

J.Phys.G G30 (2004) S51-S59. 

D. Veras, J. R. Crepp, and E. B. Ford, Formation, 

Survival, and Detectability of Planets Beyond 100 AU, 

Astrophys. J. 696 (May, 2009) 1600-1611, 

[arXiv: 0902. 2779]. 

M. Bailes, S. D. Bates, V. Bhalerao, N. D. R. Bhat, 
M. Burgay, S. Burke-Spolaor, N. DAmico, S. Johnston, 
M. J. Keith, M. Kramer, S. R. Kulkarni, L. Levin, A. G. 
Lyne, S. Milia, A. Possenti, L. Spitler, B. Stappers, and 
W. van Straten, Transformation of a Star into a Planet 
in a Millisecond Pulsar Binary, Science 333 (Sept., 2011) 
1717 , [arXiv:1108.520l[. 



